A Comparison of search templates for gravitational waves from binary inspiral 



Thibault Damour/ Bala R. lyer,^ and B.S. Sathyaprakash"^ 

^Institut des Hautes Etudes Scientifiques, 91440 Bures-sur-Yvette, France 
^Raman Research Institute, Bangalore 560 080, India 
'^School of Physics and Astronomy, Cardiff University, 5, The Parade, Cardiff, CF24 3YB, U.K. 

(Dated: February 7, 2008) 

We compare the performances of the templates defined by three different types of approaches: tra- 
- . . ditional post-Newtonian templates (Taylor-approximants), "resummed" post-Newtonian templates 

assuming the adiabatic approximation and stopping before the plunge (P-approximants), and further 

■ "resummed" post-Newtonian templates going beyond the adiabatic approximation and incorporat- 
' ing the plunge with its transition from the inspiral (Effective-one-body approximants) . The signal to 

noise ratio is significantly enhanced (mainly because of the inclusion of the plunge signal) by using 
. these new effective-one-body templates relative to the usual post-Newtonian ones for a total binary 

'~i ' mass m > 30Mq, and reaches a maximum around m ~ 80Mq. Independently of the question of the 

I plunge signal, the comparison of the various templates confirms the usefulness of using resummation 

. methods. The paper also summarizes the key elements of the construction of various templates and 

' thus can serve as a resource for those involved in writing inspiral search software. 

I PACS numbers: 04.3.0Db, 04.25.Nx, 04.80.Nn, 95.55.Ym 

> ■ 

c^ : 

O . I. INTRODUCTION 

o ■ 

■ The late stage evolution in a compact binary, when the component stars are under the influence of the 
' strong gravitational fields of each other and are moving at relativistic speeds, is dictated by the non-linear 
, dynamics of general relativity and is very difficult to model. In the early stages of adiabatic inspiral (that 
' is, when the inspiral time scale is much larger than the orbital time scale) it is possible to treat the problem 

^H. of motion perturbatively and to expand the general relativistic equations of motion and wave generation 
^ • formulae in a power-series in v/c, v being a characteristic velocity. (We henceforth use units such that 
G = c = 1.) However, the phasing of the gravitational wave (GW) signal derived from these perturbative 
• • ■ results becomes increasingly inadequate as the two bodies approach each other. The characteristic velocity 
[ Vpim) = {-KiTifpY/^ , corresponding to the peak of the detector sensitivity to the inspiral signal from a 
, binary of total mass m = mi -I- m2, is numerically equal (for initial LIGO; for which fp = 126 Hz)^ to 
?H ' Vpim) — O.125(m/Mo)^/'^. For a double neutron star system one already has Wp(2.8) = 0.176, while, for 
an archetypal {1QMq,1QMq) double black hole system one has Wp(20) = 0.340, quite close to the velocity 
corresponding to the last stable orbit (Iso): v\so ~ l/\/6 — 0.408. 

The present theoretical understanding has enabled the perturbative computation (via post-Newtonian 
expansions) of the binary orbit and GW phase to an order beyond the standard quadrupole formula. We 
shall use these u^-accurate results in this work. At present, we cannot count on the (hopeful) extension of 
the post-Newtonian (PN) perturbative calculations by another two orders, to order d^, because currently 
used PN techniques leave undetermined a physically crucial parameter entering at the level 0, IM- 
Moreover, as emphasized some time agoQ, the PN series (which is essentially a Taylor expansion in powers 
of v) is a rather poorly convergent expansion. More precisely, if one considers the PN expansion of the 
crucial GW flux (see, e.g., Fig. 3 of j^), one notices that the "y^-accurate and v^- accurate approximations 
start significantly deviating, in opposite directions, from the exact (test- mass) result when w > 0.2. As such 
relatively high values of v are typically involved in the calculation of the GW phasing {v ~ Vp{m) > 0.18 as 
soon as m > 2.8Mq), one has to worry that search templates based on a straightforward use of PN-expanded 
results might be inadequate for the detection and/or measurement of inspiral signals, especially for the more 
massive systems (m > lOAf© implies Vp{m) > 0.34) which are likely to be the first potentially detectable 
events. 

To address this crucial problem, we have been advocating [^0,0] a new philosophy for making the optimal 
use of existing PN results, namely, to use several re-summation techniques to improve the convergence of 
the PN series, before using them to compute the GW phasing. As of now, we have proposed and studied 
three successive stages in the definition and use of such re-summation techniques. First, we constructed 



It should be noted that the LIGO noise curve used in this paper is the currently best available and different from that used 
in DISl and DIS2. 
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time-domain signals, called P-approximants - starting from the standard PN Taylor representation - which 
possess better convergence properties and capture the expected analytical behaviour (poles and zeroes) of 
the relevant physical quantities quite well. We have shown that these new signal models, when compared 
with standard PN signals, are both more effectual (larger overlaps) and more faithful (smaller biases in 
the estimation of parameters) representations of some fiducial "exact" signals. Though time-domain P- 
approximants are better signal models than the standard "Taylor" approximants, they are computationally 
expensive to use in a data analysis exercise that searches for inspiral signals using hundreds of thousands 
of templates, which have to be correlated, with arbitrary time lags, with the detector output. This data- 
analysis computational cost is much reduced (thanks to the existence of efficient FFT algorithms) when one 
disposes of explicit analytic expressions for the Fourier transform of the templates. 

Second, we found j6j| explicit frequency- domain representations of P-approximants (as well as of standard 
PN templates) that are computationally inexpensive and are yet as faithful and effectual as the original 
time-domain models. This frequency-domain representation incorporates the "edge oscillations" due to the 
(assumed) abrupt shut-off of the time-domain signal occurring soon after the binary crosses the last stable 
orbit. In Q we emphasized that the signal to noise ratio of the first interferometric detectors is large enough 
for detection only for massive binary black hole systems of total mass m > 25Mq. For such systems the 
characteristic velocity corresponding to the peak of the detector sensitivity is Vp(25) ~ 0.37, which is very 
close to viso — l/'s/S, the velocity at the last stable orbit. One, therefore, expects that the first detections 
are most likely to concern massive^ systems (20Mq ^ to < AOMq) with Vp ~ viso- 

It is therefore crucial to push the re-summation techniques introduced in Q further so as to be able to 
describe not only the GW phasing during the last cycles before Iso crossing, but also during the transition'^ 
between inspiral and plunge, and during the plunge itself. Recently, Buonanno and Damour |7ij combined 
some of the re-summation techniques of Q and with a novel approach to the general relativistic dynamics 
of two-body systems 8] to devise an improved type of re-summation approach to the GW phasing of coalescing 
binaries, able to describe in more detail the transition between inspiral and plunge. This "effective one- 
body" approach is the first one that goes beyond the "adiabatic approximation", used both in standard 
(non-resummed) PN approximants, and in P-approximants. 

The data analysis of inspiral, merger and ring-down was pioneered by Flanagan and Hughes ,9J. They 
treated the problem of inspiral rather accurately but the merger/plunge was treated by assuming that about 
10% of the rest mass energy would be emitted during merger. This quite optimistic estimate was based on a 
crude model of the coalescence of maximally spinning black holes, and was arbitrarily extended to all cases. 
A similar back-of-the-envelope consideration of the ring-down amplitude let them to optimistically assume 
that about 3% of the rest mass energy would be emitted during ring-down. 

In this paper we discuss only non-spinning^ binaries and we make no ad hoc assumption about the total 
energy radiated during the merger phase. The Effective One-Body (FOB) formalism does not treat the 
inspiral and plunge phases separately. Indeed, in this formalism the plunge is seen as a natural continuation 
of the inspiral phase contributing (for equal masses) about 0.6 orbital cycles (or 1.2 GW cycles), with a 
total energy associated with the plunge around 0.7%. The energy emitted during the following (matched) 
ring-down phase is also found to be around 0.7% [Tlj . These energy losses are much smaller than the 
Flanagan-Hughes guesstimates of 10% and 3%, respectively. Consequently, it is unlikely that we will be able 
to detect the plunge phase of the EOB waveforms separately, irrespective of the mass of the system. This 
is in sharp contrast to the Flanagan-Hughes claim that the SNR contribution of the sole merger phase of 
massive black holes of total mass in the range 30-1000 Mq will dominate over the inspiral phase contribution. 
Note also that (for a source at 100 Mpc) the (merger-dominated) SNR of ^ reaches a maximum of 40 around 
TO ~ 2OOM0, while our (inspiral-dominated) SNR reaches a maximum of 8 around to ~ 8QMq. [It seems 
that most of the difference between Fig. 4 of Flanagan and Hughes and our Fig. ^ below, e.g. a factor 
3 (between their 25 and our 8) between the SNR for a 80 Mq source at 100 Mpc, comes from the huge 
difference in energy loss during merger.] 

In this paper we make a prediction that the merger phase will enhance the inspiral phase SNR by about 
10% for TO ~ 30Mq, and by about 300% for to 80Mq, but that it will not be in itself a significant signal. 
As a result our best candidate sources are still stellar mass black hole binaries of total mass in the range 30-90 



^ Note that, even for less massive systems, the necessity to capture more than 96.5% of the SNR, corresponding to a loss in the 
number of events by no more than 10%, implies that one must accurately control the phasing of the waveform at frequencies 
significantly higher than fp, corresponding to velocities significantly higher than Vp{m,). 

^ The P-approximants model this transition by a sharp cut-o ff in the signal. 

* In view of current black hole binary formation mechanisms lid , we think it likely that most of them will include only slowly 
spinning holes. 
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Mq. We also conclude that the ring-down phase is in itself not a significant contribution for m < 200Mq. 

There is one word of caution regarding the plunge signal: Even though the plunge lasts for only about 
half an orbital cycle its spectral content spreads over a large frequency range. Consequently, the number of 
frequency bins over which the signal spreads out is quite large and it is not advisable to use a non-optimal 
method to try to detect the plunge part in isolation. In fact, we believe that one of the robust predictions 
of the EOB approach (at least in the case of slowly spinning holes) is that the plunge signal is a smooth 
continuation of the inspiral one, and that one should use templates that are phase-coherent all over the 
inspiral-plus-plunge phase. We are aware of the approximate nature of the EOB results (especially beyond 
the Iso) and do not claim that the EOB waveform is the last word on the problem, but our position is the 
following: (i) in absence of comparatively accurate alternative results, it is important to study in detail the 
predictions coming from the EOB waveform, and (ii) we shall finally recommend to use a bank of filters 
which cover a large range of possibilities, with special weight being given to the best-tested "resummed" 
templates. 

We also hope that our work will give an additional incentive to numerical relativity groups toward com- 
puting waveforms which are at least as accurate (and physically complete) as the EOB one. In particular, 
let us recall that has proposed a new approach to the numerical computation of binary black hole coa- 
lescences: namely, to start the numerical evolution just after Iso crossing, i.e. at a stage where one can still 
trust the resummed PN estimate of the dynamics of two black holes, but where there is only 0.6 orbit to 
evolve before coalescence. To this aim Ref. 01 has provided explicit results for the initial dynamical data 
(positions and momenta in ADM coordinates) of this problem. However, apart from stimulating further 
thoughts on the problem (|l2l|'). we are not aware of the existence of numerical simulations implementing the 
proposal of Ref. 7], nor are we aware of other numerical work leading to explicit (non adiabatic) waveforms 
for coalescing binaries which could be compared to the EOB one. 

Data analysis groups associated with various ground-based interferometers are now finalizing the analysis 
software that will be used for GW searches in data that is expected to become available in a few years time. It 
is essential that these groups be aware of recent theoretical developments and of their respective merits so as 
to take the best advantage of the current knowledge in writing their software. With this view in mind the aim 
of this paper is two- fold: First, we wish to compare the performances of the templates defined by the three 
different types of approaches mentioned above (traditional "non-resummed" PN templates, "resummed" PN 
templates assuming the adiabatic approximation and stopping before the plunge, and further "resummed" 
PN templates going beyond the adiabatic approximation and incorporating the plunge with its transition 
from the inspiral.) Second, in view of the fact that the original publications are quite complex and 

technical, we wish to summarize in a more accessible manner the key elements of the techniques introduced 
there (and re-used, with further inputs, in 0-) The present work should serve as an easily accessible resource 
for data analysis groups. Readers interested in a more detailed understanding of our general approach are 
referred to 5, 6, 7, 8] for motivation, formalism, logical reasoning, exhaustive tests and further discussion of 
the new signal models. 

II. TIME-DOMAIN PHASING FORMULAS IN THE ADIABATIC APPROXIMATION 

In searching for GW from an inspiralling compact binary we are faced with the following data analysis 
problem: On the one hand, we have some (unknown) exact gravitational waveform h^{t]\k) where Afc, 
k = 1, . . . ,n\, are the parameters of the signal (e.g., the masses mi and TO2 of the members of the emitting 
binary). On the other hand, we have theoretical calculations of the motion of 0,|2>S^3) and gravitational 
radiation from [T^ ITsj binary systems consisting of neutron stars (NS) or black holes (BH) giving the 
PN expansions of an energy function E{x = w^), which is related to the total relativistic energy Etot via 
^-tot = i'Tti +m2){l + E), and a GW luminosity (or "flux") function T{v). Here, the dimensionless argument 
V = IS an invariantly defined "velocity" related to the instantaneous GW frequency F (= twice the 
orbital frequency) hy v = {TrmF) a . Given PN expansions of the motion of and gravitational radiation from 
a binary system, one needs to compute the "phasing formula", i.e. an accurate mathematical model for the 
evolution of the GW phase^ = p[t; A;], involving the set of parameters {A^} carrying information about 
the emitting binary system. In the adiabatic approximation the phasing formula is easily derived from the 
energy and flux functions. Indeed, the standard energy-balance equation dEtot/dt = —T gives the following 



^ We work within the "restricted" waveform approximation which keeps only the leading harmonic in the GW signal. 
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parametric representation of the phasing formula: 



°' , E'{v) , „ , 3 E'{v) 



t{v)^U.^f + m dv—^, = 0,.cf + 2 / dw-"—^, (2.1) 

where t-^ci and (/)rof are integration constants and v-[ci an arbitrary reference velocity. [It is sometimes 
convenient, though by no means necessary, to take as Uj-cf the velocity v at the last stable orbit (see below).] 
^From the view point of computational purposes it is more efficient to work with the following pair of coupled, 
non-linear, ordinary differential equations (ODE's) that are equivalent to the above parametric formulas: 

d(l> 2t;3 _ ^ dv ^ T{v) _ ^ ^2 2^ 

dt m ^ dt mE'{v) 

We shall see later that, for massive systems, the adiabatic approximation fails and one must replace the 
two ODE's by a more complicated ODE system. We now turn to a discussion of what is known about the 
two functions E{v) and T{v) entering the phasing formula and how that knowledge can be improved. 

A. T-approximants 

We denote by Et„ and Tt„ the n'^-order^ "Taylor" approximants (as defined by the PN expansion) of 
the energy and flux functions: 



EtM = ENix)J2Ek{ri)x\ Tt^{x) = Tn{x) 



k=0 



^.Ffc(r,)«'=+^Lfe(77)log(z;/z;o)z;'= 



.k=0 k=6 



(2.3) 



1 32 

where, En{x) — — rjx, Tn{x) — — ifx^. (2-4) 

2 5 

Here the subscript N denotes the "Newtonian value", rj = raim2/rn? is the symmetric mass ratio, and wq 
is a fiducial constant to be chosen below. In the test mass limit, i.e. 77 — > 0, E{x) is known exactly, from 
which the Taylor expansion of E t„ (f , ) in Eq. (|2.3|l . can be computed to all orders. In the 77 liinit, the 
exact flux is known numerically [la, [l3 and the Taylor expansion of flux in Eq. (|2.3I) is known (13, lla, 113 
up to order n = 11. On the other hand, in the physically relevant case where 77 is finite, the above Taylor 
approximants are known |l4l ITsI l only up to five-halves PN order, i.e. n = 5. Recently the energy function 
has been computed up to 3PN order i.e. rt = 6, though with the presence of an unknown parameter |^. 
All the completely known coefficients in the expansions are enlisted in Tabled 

The problem is to construct a sequence of approximate waveforms h^{t\ Xk), starting from the PN expan- 
sions of E(v) and J-'{v). In formal terms, any such construction defines a map from the set of the Taylor 
coefficients of E and J- into the (functional) space of waveforms. Up to now, the literature has considered 
(one of) the map(s), say T, 

iET^.TrJ^hlit.Xk), (2.5) 

obtained by inserting the successive Taylor approximants into the phasing formula 0, . For brevity, we 
often refer to these "Taylor" approximants as " T-approximants". It must also be emphasized that even 
within this Taylor family of templates, there are at least three ways of proceeding further, leading to the 
following three inequivalent constructs: 

tl.) One can retain the rational polynomial ^t„/Et„ as it appears in Eq. (|2.2|) and integrate the two ODE's 
numerically. We shall denote the phasing formula so obtained as cf)^^ (t) : 

^=0, ^ + 4-^=0. (2.6) 



dt m ^ dt mE'j. (y) 



The label n always refer to an approximant accurate up to v" = a;("/2) included. 
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t2.) One can re-expand the rational function Tt^/Et„ appearing in the phasing formula and truncate it at 
order w", in which case the integrals in Eq. H2.1() can be worked out analytically, to obtain a parametric 
representation of the phasing formula in terms of polynomial expressions in the auxiliary variable v 

n n 

'/'tJ(«) = + E = 41 + tl{v) E (2-7) 

fc=0 fe=0 

where a superscript on the coefficients (eg. indicates that v is the expansion parameter [as is 
explicit from Table ITU the coefhcient of (/)^ include in some cases, a logu dependence]. 

t3.) Finally, the second of the polynomials in Eq. 12.7|l can be inverted and the resulting polynomial for v 
in terms of t can be substituted in [v) to arrive at an explicit time-domain phasing formula 

n n 

-^rjW - '/'^cf +'^^E'^fe^'' = Fl,Y.Fl6\ (2.8) 

where 9 — [ri{t-^cf ^ <)/(5to)]^^/® and F = d(l>/2ndt = v'^/{TTm) is the instantaneous GW frequency. 
The coefficients in these expansions are all listed in Table Ull 



B. P-approximants 

Before defining new "resummed" energy and flux functions with improved performances we digress for a 
brief reminder of Pade re-summation, which is a standard mathematical technique used to accelerate the 
convergence of poorly converging power series. Let Sn{v) = ao + aiv + ■ ■ • + anv"" he a truncated Taylor 
series. A Pade approximant of the function whose Taylor approximant to order is Sn is defined by two 
integers m, k such that m + k — n. If T„[- • • ] denotes the operation of expanding a function in Taylor series 
and truncating it to accuracy (included), the P™ Pade approximant of Sn is defined by 

^"W^^TT; T4Pi{v)]^Sn{v), (2.9) 
Dk[vj 

where Nm and Dk are polynomials in v of order m and k respectively. If one assumes that Dk{v) is normalised 
so that Dk{Q) — 1; i.e. Dk{v) — 1 + (/i u + • • • , one shows that Pade approximants are uniquely defined by 
Eq. (|2.9(l . In many cases the most useful^ Pade approximants are the ones near the "diagonal", m = fc, i.e. 
P™ \in — 2m is even, and P™'*'^ or P™_|_i if n = 2m + 1 is odd. In this work we shall use the diagonal (P™) 
and the "sub-diagonal" (P™_|_]^) approximants which can be conveniently® written in a continued fraction 
form '20'|. For example, given 5*2 (w) = + ai v -\- 02 one looks for 

pi.. Co Co(l +C2f) , . 

The continued fraction Pade coefficients of a series containing six terms, i.e. S^{v), are given by 

ai a2 , ai ao(aia3 - a\) 
Co = oo, ci = - — , C2 = -— -I- — , C3 = ^ 2 



C4 



ao' ai flo' ai (0^ — 0280) 

CoCl(c2 -I- CiY -I- CoCiC2C3(c3 + 2c2 + 2ci) - 04 

C0C1C2C3 

((C2 +Cl)2 +C2C3)2 (C4+C3+C2+Ci)2 flg 



C5 = • (2-11) 

C2C3C4 C4 C0C1C2C3C4 

In and we introduced several techniques for "re-summing" the Taylor expansions (in powers of v) 
of the energy and flux functions. Starting from the PN expansions of E and T , in DISl we proposed a 



The rare theorems dealing with the Pade technique concern the convergence of "near-diagonal" Pade approximants, i.e. 
m — > 00 with \m — k\ fixed. 

A convenience of this form is that the n-th continued-fraction coefficient c„ (see below) depends only on the knowledge of 
the PN coefficients of order < n. 
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new class of waveforms, called P-approximants, based on two essential ingredients: (i) the introduction, on 
theoretical grounds, of two new, supposedly more basic and hopefully better behaved, energy-type and flux- 
type functions, say e{v) and /(w), and (ii) the systematic use of Pade approximants (instead of straightforward 
Taylor expansions) when constructing successive approximants of the intermediate functions e(f), f{v). 
Schematically, our procedure is based on the following map, say "P" : 

{Et^.TtJ ^ (eT„,/Tj (ep„,/p„) ^ (i?[epj, ^[ep„ , /pJ) ^ hP{t,Xk). (2.12) 
Our new energy function e{x), where x = v^, is constructed out of the total relativistic energy Stotl^;) using 



e{x) 



2rniTO2 



- 1. 



(2.13) 



The function E{x) entering the phasing formulas is the total energy per unit mass after subtracting out the 
rest mass energy: E{x) = [E^otix) — m]/m and is given in terms of e{x) by 



E{x) = [l + 2r; ( + ei^) - l) 



1/2 



^ dE 

' dx 



r]e'{x) 



2[1 + E{x)] .yi + e{x) 



(2.14) 



Note that the quantity E'{v), needed in the phasing formula, is related to dE{x)/dx via E'{v) = 2vdE{x)/dx. 
In the test-mass limit e(x) and dE(x)/dx are known exactly: 



.0(2^) = -X 



I -Ax 
1 - 3a;' 



E''='\X)=TJ 



1 - 2x 



- 1 



dx 



■q (1 — 6a;) 
"2 (1 -3a;)3/2- 



(2.15) 



The rationale for using e{x) as the basic quantity rather than Etot{x) can be found in 5|; here we note the 
following two points: (1) In the test mass case e{x) is meromorphic in the complex a;-plane, with a simple 
pole singularity, while the function E{x) is non-meromorphic and exhibits a branch cut. (2) Secondly, in 
the test mass case, the Pade approximant of eT2„(2;), for n > 2, yields the known exact expression including 
the location of the Iso and the pole. Therefore, the function e(x) is more suitable in analyzing the analytic 
structure than is E{x). In the comparable mass case, under the assumption of structural stability between the 
case 77 — > and the case of finite 77, one can expect the exact function e{x) to admit a simple pole singularity 



,) . We do not know the location of this singularity, but Pade approximants are 



on the real axis oc (a; — Xpoi, 

excellent tools for giving accurate representations of functions having such pole singularities 

Our proposal is the following: Given some usual Taylor approximant to the normal energy function, 
Et^„ = -\rix {I + Eix + E2X^ + ■ ■ ■ + E^ a;"), one first computes the corresponding Taylor approximant 
for the e function, say 



-X > CkX 

k=0 



(Remember that we consistently label the successive approximants by the order in velocity; e.g. 
accurate object has the label 4.) Then, one defines the Pade approximant of eT2„(a;)^ 



eP2, 



E 

Lfe=0 



CkX 



(2.16) 



a 2PN- 



(2.17) 



where e = or 1 depending on whether n = 2m + e is even or odd. We shall call the continued fraction 
Pade coefficients of ep^^ as ci, C2, • • • , (Note that co = 1). They are given in terms of efe in Table HI Given a 
continued fraction approximant ep2„(a;) of the truncated Taylor series eT2„ of the energy function e[x) the 
corresponding Ep^^(x) and dEp^^(x)/dx functions are obtained using formulas in Eq. (|2.14|) by replacing 
e{x) on the right hand side with ep2^{x). 

Apart from using it to improve the convergence of the PN series, Ref. ^ has also proposed to use the 
Pade-resummed function ep2„ (a;) to determine the location of the Iso, the Pade estimates of the Iso being 
defined by considering the minima of ep2„(a;). In contrast, in the Taylor case one must, for consistency. 



More nreciselv. p, d_ (x) is —x times the Pade annroximant of —x ^e.rr^- (x). 
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use the minima of Et„ (v) to define the locations of the Iso. We have confirmed that in the test mass case 
this Pade-based method yields the exact result at orders and beyond while the corresponding Taylor- 
based method [considering the minima of Et„{v)] gives unacceptably high estimates of viso, i-e. of the 
GW frequency at the Iso . In the finite 77 case, the Pade-resummed predictions are in good qualitative, 
(and reasonable quantitative) agreement with the more recent predictions based on the "effective-one-body" 
approach!^. The location of the Iso's for the various approximations together with the location of the light 
ring [i.e. the pole singularity in ep2^{x)] are also tabulated in Tabled 

Having defined a new energy function, we move on to introduce a new flux function. The aim is to define 
an analytic continuation of the flux function to control its analytic structure as also to handle the logarithmic 
terms that appear in the flux function in Eq. H2.3|) . Factoring out the logarithmic terms is what allows us 
to use standard Pade tec hniq ues effectively in this problem. 

It has been pointed out [13 that the flux function in the test mass case J^{v; rj — 0) has a simple pole at the 
light ring = 1/3. It has been argued that the origin of this pole is quite general [cf. 5], discussion following 
Eq. (4.3)] and that even in the case of comparable masses we should expect to have a pole singularity in J^. 
However, as already pointed out, the light ring orbit in the i] ^ case corresponds to a simple pole a;poic('7) iu 
the new energy function e{x; r/). Let us define the corresponding (invariant) "velocity" Wpoic('7) = \/ Xpoicii])- 
This motivates the introduction of the following "factored" fiux function, f{v]ri) 

f{v;v) = (l-«/i.poic) Hv;v)- (2.18) 

where J-'{v) = !F{v) / Tn{v) — 5J-'{v)/ {32ri^v^'~'), is the Newton-normalised fiux. Note that multiplying by 
1 — u/fpoic rather than 1 — (u/wpoie)^ has the advantage of regularizing the structure of the Taylor series of 
f{v) in introducing a term linear in v, which is absent in the fiux function in Eq. H2.3|l (cf. Table Hlj) . Three 
further inputs will allow us to construct better converging approximants to f{v). First, it is clear (if we think 
of V as having the dimension of a velocity) that one should normalize the velocity v entering the logarithms 
in the fiux function in Eq. (|2.3|) to some relevant velocity scale vq. In the absence of further information the 
choice Vq = visoiv) seems justified (the other basic choice vq — Upoic is numerically less desirable as v will 
never exceed viso and we wish to minimize the effect of the logarithmic terms). A second idea, to reduce the 
problem to a series amenable to Padeing, is to factorize the logarithms. This is accomplished by writing the 
/ function in the form 

/tJ«;?7)= (l + E^^^'^ln^J (2.19) 

\ k—6 / \k—0 / 

where coefficients /j. are /q = 1, fk+i = ^k+i—^k/v^°^^ and 1^ are constants determined from the coefficients 
of by relations (analogous to) Eq. (4.9) of DISl^*^. Finally, we define Pade approximants to the factored 
fiux function f{v) as 

V k=6 ^lso('/)/ Vfc^O / 

where v^^iv) denotes the Iso velocity (= y/xiso) for the t;"-accurate Pade approximant of e(x), and where 
Pm+e denotes as before a diagonal or sub-diagonal Pade with n = 2m + e, e = or 1. The corresponding 
approximant of the fiux J-{v) is then defined as 

Tpjv;rj)= (1-^^] fpAv;v), (2.21) 

where f^oieC'?) denotes the pole velocity defined by the w"-Pade of e{x). In the test mass case the exact 
location of the pole and the Iso are Xpoio = 1/3 and x\so = 1/6, respectively (cf. Table CJ. We shall denote 
the continued fraction Pade coefficients of fp^{v) by dk- They can be found in terms of fk using Eqs. H2.11|l . 
At present, the most accurate estimate of the flux would be ^Pn, defined by using the known ?7-dependent 
coefficients /fe for /c < 5, and the test- mass values of fk and Ik for /c > 6. 



The variables 1^ and /j. used here are equal to the variables 1^ and /j. used in DISl. They are 'hatted' here as a reminder 
that they represent coefficients of Newtonian-normalised quantities. The coefficients Aj. and B/. appearing in the definition 
of lu are rommited in R.ef. 
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III. FREQUENCY-DOMAIN PHASING FORMULAS IN THE ADIABATIC APPROXIMATION 



The time-domain P-approximant waveforms discussed above are computationally intensive to use in a 
full-scale data-analysis. Recently, in DIS2 we have constructed frequency-domain representations of the 
P-approximants which are 10-50 times faster to compute than their time-domain analogues but are yet as 
accurate. This increases the usefulness of P-approximants in data analysis. 

The Fourier representations are based primarily on a newly derived improved version of the stationary 
phase approximation appropriate to time-truncated signals. P-approximants cannot be modeled using the 
standard stationary phase approximation over the entire frequency domain. Indeed, close to the last stable 
orbit, where the inspiral phase terminates, one requires a modification of the stationary phase approximation. 
In simple terms, Ref.Q found a way of taking into account the effect of an assumed abrupt termination of 
the waveform near the last stable circular orbit by introducing simple modifications to the usual stationary 
phase approximation. We present only the final results here; the interested reader is referred to DIS2 for 
details. Note that the results summarized below are quite general and can be applied to a generic chirp 
signal which shuts off abruptly (i.e. on a time scale ^ F^^). 

We begin with a discussion of the usual stationary phase approximation for chirp signals. Consider a 
signal of the form, 

/i(t)==2a(i)cos0(t)=a(t)[e-*'^(*)-fe*'^(*)j , (3.1) 

where (p{t) is the implicit solution of one of the phasing formulas in Eq. I|2.6|l . Eq. (|2.7I) or Eq. (|2.8|) for some 
choice of functions E' and !F . 

The quantity 2-KF{t) — d(j){t)/dt defines the instantaneous GW frequency F{t), and is assumed to be 
continuously increasing. (We assume F{t) > 0.) Now the Fourier transform h{f) of h{t) is defined as 

/"Oc poo 

h{f)= / dte^'^'f'hit) ^ / dta{t) e^^^ft-m ^^^^ft+m (3 2) 



The above transform can be computed in the stationary phase approximation (SPA). For positive frequencies 
only the first term on the right contributes and yields the following usual SPA: 

^uspa^j) ^ a{tf) ^»[^j(t,.)-V4] ^^.(^^ ^ 27r/i- 0(i), (3.3) 



and i/ is the saddle point defined by solving for dipf{t)/dt = 0, i.e. the time tf when the GW frequency 
F{t) becomes equal to the Fourier variable /. In the (adiabatic) approximation where Eqs. (|2.1() hold, the 
value of tf is given by the following integral: 

tf^t,,f + m -TEF^dv, (3.4) 
where Vf = (ttto/)^/'^. Using tf from the above equation and (j){tf) in Eq. (|2.1|) one finds that 



4,f{tf) = 27r/W - (/-ref + 2 J {vj- 



"rot 



E'{v) 



dv. (3.5) 



The big computational advantage of Eq. (|3.5|) [with respect to its time-domain counterpart, Eq. H2.1|l ]. is 
that, in the frequency domain, there are no equations to solve iteratively; the Fourier amplitudes are given 
as explicit functions of frequency. 

In the Fourier domain too there are many inequivalent ways in which the phasing ^f can be worked out. 
Here we mention only the most popular: 

fl.) Substitute (without doing any re-expansion or re-summation) for the energy and flux functions their 
PN expansions or the P-approximants of energy and flux functions and solve the integral in Eq. H3.5() 
numerically to obtain the T-approximant SPA or P-approximant SPA, respectively. 

f2.) Use PN expansions of energy and flux but re-expand the ratio E'{v)/T{v) in Eq. 13.51) in which case 
the integral can be solved explicitly. This leads to the following explicit, Taylor-like, Fourier domain 
phasing formula: 

5 

Xj;f{tf) = 2^/trcf - <^rcf +TAr^ffc(7rm/)('=-^)/3 (3.6) 

fc=0 

where tu are the chirn narameters listed in Table ITTl 
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Eq. is one of the standardly used frequency-domain phasing formulas. Therefore, we shall use that 

as one of the models in our comparison of different inspiral model waveforms. We refer to it as "type-f2" 
frequency-domain phasing. 

Just as in the time-domain, the frequency-domain phasing is most efficiently computed by a pair of coupled, 
non-linear, ODE's: 

rather than by numerically computing the integral in Eqs. I|3.5|l . 

Next we correct the performance of the usual SPA by including edge corrections arising as a consequence 
of modeling the time-domain signal as being abruptly terminated at a time t = tmax (time-truncated chirp) 
when the GW frequency reaches F = Fmax- In practice, we expect that Fmax will be of order the GW 
frequency at the Iso. However, we prefer to leave unspecified the exact value of Fmax- The idea is to use 
-Fmax as a free model parameter, to be varied so as to maximize the overlap (see Section 01 between the 
^max-truncated template and the real signal. Such time-truncated signals can be represented as: 

h{t) = 2a{t) cos0(t)6'(t,„ax - i), (3.8) 

where 9 denotes the Heaviside step function, i.e. 0{x) = if x < and 6{x) = 1 when x > 0. The effect 
of this time-windowing has been modeled in DIS2 and the result is that the Fourier transform of such a 
time-truncated signal can be accurately represented in the two regions / < i^max and / > i^max, by 



/<i^.„ax : /i'r(/)=C(C<(/))^i^e^['^/(*/)--/4l, 



f>F,^,^ : hTif)=C{C>{.f)) "^^"^'^^ exp 



^i^(imax) 



*^/'/(tmax) + ^ ^^{^,,^'°7'' ~ ^'^Z'* 



(3.9) 



where the label 'ispa' stands for improved SPA, C is essentially the complementary error function, 2C((^) = 
erfc (e''^/'*C) , and (J is computed using 

/ < Fmax : C< = ~[^f{tf) - ^/(tmax)]'/', 

/ > Fmax : C>(/) = ^^^-^---\ (3.10) 



max } 

The error function needed in calculating C(/) may be numerically computed using the NAG |22 library 



S15DDF. Note the denominator \j F(tmax) entering Eqs. (|3.9|l and (|3.1U|I . We define a generic time-truncated 
signal as a chirp for which this denominator is finite (neither infinite, nor zero). [These signals were called 
"Newtonian- like" in Ref. 6]. We do not keep this name here to emphasize that the results Eqs. H3.9|l and 
(|3.1UI) apply also to relativistic models, as long as F(tmax) is finite.] 



The exceptional (non-generic) case where y F(imax) becomes infinite arises if one tries to keep using the 
simple adiabatic phasing approximation up to the last stable orbit defined by the corresponding (approxi- 
mate) energy function E(v). This exceptional case can also be dealt with at the price of a more complicated 
modification of the usual stationary phase result (see Ref. ^ for details). We do not enter into details here 
because the recent work 7] on the transition between the adiabatic inspiral and the plunge has shown that 
the adiabatic approximation breaks down just before the Iso, and that F(t) never becomes infinite. 



IV. TIME-DOMAIN PHASING FORMULA BEYOND THE ADIABATIC APPROXIMATION 

In the above, we restricted ourselves to the standard "adiabatic approximation" , where one estimates 
the phasing by combining the energy-balance equation dEtot/dt = —T with some resummed estimates for 
the enersi and flux as functions of the instantaneous circular orbital frequency. Recently, Buonnano and 
Damour |3 have introduced a new approach to GW from coalescing binaries which is no longer limited to 
the adiabatic approximation, and which is expected to describe rather accurately the transition between 
the in.snira.l a.nd the nlunffe. and to erive also a, fir.st estimate of the following nlunffe siffna.l. The a.nnroar.h 
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of is essentially, like 0, @i a re-summation technique which consists of two main ingredients: (i) the 
"conservative" (damping-free) part of the dynamics (effectively equivalent to the specification of the E{v) 
in the previous approaches) is resummed by a new technique which replaces the two-body dynamics by an 
effective one-body dynamics Q , and (ii) the "damping" part of the dynamics [equivalent to the specification 
of the ^{v)] is constructed by borrowing the re-summation technique introduced in In practical terms, 
the time-domain waveform is obtained as the following function of the reduced time t = t/m: 



h{i)=Cvl{i)cos{(l)Gw{i)). 



dip 
di 



The orbital phase (p{i) entering Eq. (|4.1|) is given by integrating a system of ODE's: 

dr dH , 



dt 

dip 
di ' 

dpr 

di 

dp^ 

di 



(4.1) 

(4.2) 

(4.3) 

(4.4) 
(4.5) 



The reduced Hamiltonian H (of the associated one-body problem) is given, at the 2PN approximation, by 



H{r,pr,p^) = - A 



1 + 27] 



2 277 , , 1 

where A{r) = 1 - - + , B{r) = —— 
r r'^ A(r) 



1 _ ^ 



The 3PN version of H has been recently obtained The damping force J-^p is approximated by 
where Tp^{vi^) — ^ rj^ ^p„{vlj) is the flux function used in P-approximants discussed above. 



(4.6) 
(4.7) 

(4.8) 



The system Eq. (jIVII allows one to describe the smooth transition which takes place between the inspiral 
and the plunge (while the system H2.2|) becomes spuriously singular at the Iso, where E'{v\so) = 0). Ref.0 
advocated to continue using Eqs. (|1V|) after the transition, to describe the waveform emitted during the plunge 
and to match around the "light ring" to a "merger" waveform, described, in the first approximation, by the 
ringing of the least-damped quasi-normal mode of a Kerr black hole, [see Eq. (6.2) of This technique 
is the most complete which is available at present. It includes (in the best available approximation and for 
non-spinning black holes) most of the correct physics of the problem, and leads to a specific prediction for the 
complete waveform (inspiral + plunge + merger) emitted by coalescing binaries. Because of its completeness, 
we shall use it as our "fiducial exact" waveform in our comparison between different search templates. 

The initial data needed in computing this effective one-body waveform are as follows: In gravitational wave 
data analysis we are normally given an initial frequency /o (wq = 7i'™/o) corresponding to the lower cutoff of 
a detector's sensitivity window, at which to begin the waveform. The initial phase of the signal will not be 
known in advance but in order to gauge the optimal performance of our approximate templates we maximise 
the overlap (see Section|Vll over the initial phases of both the fiducial exact signal (i.e., the effective-one-body 
waveform) and the approximate template. The general analytical result of this maximisation was discussed 
in Appendix B of DlSl. In the terminology of this Appendix, these fully phase-maximised overlaps were 
called the best overlaps (they are given by Eq. (B.ll) of DISl). As discussed in Appendix B of DISl, there 
are two distinct measures of the closeness of two signals: the best overlap (maximised over the phases of 
both the template and the exact signal), and the minimax overlap (maximised over the template phase, with 
the worst possible exact phase). In an investigation such as ours (where we are interested in the optimal 
mathematical closeness between different signals), the best overlap is the mathematically cleanest measure 
of closeness of two families of templates, and we shall use it here. In addition, we shall also maximise over 
the other temnlate narameters fin narticular. the masses! to set an intrinsic measure of the closeness of 
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two families of templates. Note that the resulting fully maximised overlaps are different from the maximised 
ambiguity function of Ref . l23l| and the fitting factor of Ref. |24| . The latter (identical, but given different 
names by different authors) quantities are well-defined measures of the closeness of two signals only within 
the (simplifying) approximation where signals in quadrature are orthogonal. This is, however, not the case 
for the signals we consider, and at the accuracy at which we are working. 

For the computation of the best overlaps, it is sufficient to construct two signal waveforms, and two 
template waveforms, one with phase equal to and another with phase equal to tt/2. The rest of the initial 
data (ro,p°,p°) are found using 



'0 



1 + 2?7(V^-1) 



1 







UJ, 







0, £ 



377 



3rl 



5?7 



1/2 



Ciro,pl){dpl/dro) 



where z(r) and C{r,pip) are given by 



z{r) 



3r^ + 5ri 



A^{r) 



(4.9) 



(4.10) 



The plunge waveform is terminated when the radial coordinate attains the value at the light ring r\j- given 
by the solution to the equation. 



The subsequent "merger" waveform is constructed as in Ref. Q • 



(4.11) 



V. RESULTS AND CONCLUSIONS 

In this section we compare the performances of various signal models by choosing as fiducial "exact" 
signal model the effective one-body waveforms discussed in the previous Section. An important yardstick 
for comparing different waveforms is the overlap: Given two waveforms h and g their overlap is defined as 



0{h,g) = 



{h,g) 



{h,h)'/' {g,g) 



1/2- 



In the above equation the scalar product (, ) is defined as 



(Kg) 



df 



Shif) 



Hf)9*if) + C.C. 



(5.1) 



(5.2) 



where C.C. denotes complex conjugation and Sh{f) is the one-sided detector noise spectral density 



one — sided 



= 25: 



two— sided 



leading to the factor 2 in Eq. H5.2|l . compared to the definition used in 0, where we 



always use the two-sided noise). See the Appendix below for the noise performances of the various detectors. 

Firstly, in Fig. ^ we compare the signal-to-noise ratios (SNRs), expected in GEO, LIGO and VIRGO, for 
equal mass binaries located at 100 Mpc, when detecting an "exact" signal h by means of a bank of templates 
k: 



N 



(5.3) 



Thick lines plot the SNR obtained when k = h, i.e. when the template perfectly matches our fiducial exact 
waveform (i.e. the effective one-body waveform including its "ringing tail"), and thin lines show how that 
gets degraded when we use for k the best post-Newtonian template T-''^ [cf. Eq. (j3.6|l ] truncated at the test- 
mass Iso, = 44OOM0/TO Hz, assuming still that the true signal h is the one-body effective waveform. 
(As usual, see, e.g.. Section IV A of 0, we averaged over all the angles.) The overlaps 0{k, h) are maximised 
over the time lag and the two phases (as explained in the previous section), as well as over the two individual 
masses mi and 7712.^^. The greater SNR achieved by effective one-body waveforms for higher masses, as 



The plots are jagged because we have computed the SNR numerically by first generating the waveform in the time-domain 
and then iLsine its discrete Fourier transform in Eo.l5.l-il 
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compared to Fig 1 of DIS2, is due to the plunge phase present in these waveforms. We have checked that 
the final merger signal, modelled as a quasi-normal-mode, has a numerically insignificant effect in both the 
SNR and the overlaps: the overlap between our fiducial exact waveform and an effective one-body waveform 
minus quasi-normal modes is greater than 0.98 for double black holes of masses smaller than (4O,4O)M0. 

The origin of the enhancement in SNR is easily understood. The plunge waveform begins around / ~ i^iso 
and lasts until / ~ 2Fiso . The mass of a binary whose last stable orbit velocity is equal to the characteristic 
velocity Vp of LIGO's peak sensitivity is m = 34.8Mq. Therefore, effective one-body waveforms from binaries 
of masses in the range 35 Mq < m < 70Mq have larger SNRs for LIGO, than the usual SPA models. It 
is, therefore, crucial to go beyond the waveforms given in the adiabatic approximation to take advantage of 
these higher SNRs for larger masses. Indeed, the SNRs being as high as 8, for one detector and a source at 
100 Mpc, a network of four detectors (two LIGOs, GEO and VIRGO) will be able to reliably search for such 
systems as far as 150-200 Mpc. 

Next, in Table ITTll we show the fully maximised overlaps of effective one-body waveforms with signal models 
T" (column 2), r*^ (column 3), T^^ (column 4), T^^ (col umn 5) and P (column 6) for four typical binaries. 
The time-domain T*'^-approximants are terminated when F = and the other approximants are terminated 
when F{t) — fiso, the Iso being determined consistently using E'j^{v) = 0, where Ea{v) is the corresponding 
approximate energy function. The overlaps are computed with the (initial) LIGO noise curve given in the 
appendix. 

In addition to maximization over the time lags and the phases as explained earlier, we also maximize over 
the masses of the approximate waveforms keeping the masses of the exact waveforms fixed. We note that 
none of the models have good overlaps with the "exact" one for heavier mass binaries. This is as expected 
since it is for heavier masses that the characteristic plunge phase makes a significant difference between 
the approximate and the exact models. The relative performances of the 2PN and 2.5PN Taylor-templates 
depends on the choice of the scheme used as is evident from columns three and four in the Table. This is 
consistent with the results of Ref. psj . (The numerical results must not be compared since they quote values 
for the "advanced" LIGO.) Some Taylor models are "effectual" (large maximized overlaps), but at the cost of 
high biases in the parameters (i.e, in the terminology of they are not "faithful"). For example in the case 
of a fiducial exact (1.4,10) Mq system T^^ reaches 0.9452 for mass values (0.8,20) Mq! Thus T-approximants 
in the time-domain are significantly inferior to the other models, both for their erratic convergence properties 
( the w^-accurate templates being worse than the w^-accurate ones) and for their poor parameter estimation 
performance. 

To further explore the performance of the various models we plot in Fig. [21 their intrinsic frequency 
evolutions in the LIGO band, i.e. we plot F/F"^ versus F. The plot corresponds, in the fiducial exact 
case, to a binary black hole of (10,10) Mq. For the approximate models we exhibit the frequency evolution 
of the system that achieves the maximum overlap. As expected, the maximum overlap is obtained for 
template parameters such that the intrinsic frequency evolution of the template waveform is "tangent" to 
the exact one, near the maximum sensitivity of the detector. This can always be achieved by fitting the mass 
parameters. The question is whether such a local "tangency" ensures a sufficiently good "global" agreement. 
For instance, we note that the T*^''^ 2.5 PN models fare poorly in globally mimicking the frequency evolution 
of the exact waveform. This is consistent with their returning the worst overlaps of all. On the other hand, 
even though the T^^ models do not reproduce the exact model over as large a range as the P-approximants, 
they achieve nearly as large overlaps as the P-approximants, because they can be made (by optimizing 
the masses) to agree well with the exact model over most of the sensitive part of the LIGO band. The P- 
approximants are able to mimic the "exact" evolution the best with little bias in the masses but, being based 
on the adiabatic approximation, they fail to capture the smooth transition to plunge-'^^. The filters using 
the effective one-body approach go beyond the adiabatic approximation and include a smooth transition 
to plunge and merger. They, therefore, supersede the adiabatic-limited P-approximants. This difference 
between the two re-summed versions of binary signal models is important for masses larger than about 20 

Note that, following we have generated the effective-one-body model using the adimensional time 
i = t/m. [This trivially extends beyond the 2PN approximation Q.] It has been recently emphasized 
|26j| (in the context of the T- and P-approximants, where one can also simplify formulae by working with 
t = t/m) that there are many computational advantages in working with such adimensionalized time models. 
Indeed, the phase evolution becomes completely independent of the total mass of the system. This, together 



We have confirmed that the P-approximants return the best overlaps when truncated at the P— defined Iso. Maximizing over 
a cutofl^ frequency smaller than the P— defined Fi^^ (which turns out to be higher than the efTective-one-body-defined one) 
does not improve the overlaps. 
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with the fact that the evolution can be computed from a system of ODE's, makes the computation and 
storage of templates required in a search for binary black holes and neutron stars in interferometer data 
computationally inexpensive as compared to the conventional method that uses a 2-dimensional lattice of 
templates. 

We are currently estimating the effects of unknown parameters in the 3PN motion and wave generation 
■ The extension of the type of work presented here to go beyond the restricted post-Newtonian approach, 
and also to include the effects due to spin and eccentricity, needs to be systematically investigated. 

To conclude: we believe that many of the new technical tools developed in 0| and briefly summarized 
above are useful ingredients for constructing effectual and fast-computed inspiral templates. For example: 
(i) the specific Pade-based resummation of the GW flux introduced in ^ is an important ingredient of the 
construction in (3] of an accurate non-adiabatic waveform and, (ii) the improved SPA technique derived in 
could be used to derive analytical approximations to the frequency-domain version of these effective-one- 
body waveforms. In view of our ignorance of the "exact" waveform emitted near and after the Iso crossing, 
the best strategy is probably to construct a bank of templates which cover a large range of possibilities with 
special weight being given to the templates incorporating the best tested re-summation methods (such as 
P-approximants, and the effective-one-body-approach0). Because of the admittedly quantitatively rough, 
but plausibly qualitatively correct, description of the plunge signal given by the EOB approach, we also 
recommend to include some sort of multi-parameter template which qualitatively looks like Fig. 12 of 0, 
but which introduces some flexibility both in the phasing evolution during the plunge, and in the location 
of the matching to the ring-down (with the possible inclusion of several quasi-normal modes). Finally, we 
emphasize the importance of modelling the transition to the plunge, and of including the signal emitted 
during the plunge: this leads to a very significant enhancement of the signal-to-noise ratio, from about 4.5 
to 8, for a source at 100 Mpc. 
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APPENDIX A: NOISE POWER SPECTRUM OF INITIAL INTERFEROMETERS 

In this short appendix we list the expected one-sided noise power spectral densities of the various ground- 
based interferometers (Table llVjl and plot the effective noise hn = \/ fShif) in Fig. |21 
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TABLE I: Taylor coefficients of the energy functions Et„{x) and eT„(s) and the corresponding location of the Iso and 
pole. As there are no terms of order we have exceptionally chosen (for this Table only) the expansion parameter 

to be a; = rather than v. In all cases the k — coefRcient is equal to 1, the last stable orbit is defined only for 
A: > 1 in the case of T-approximants and for fe > 2 in the case of P-approximants and N denotes the "Newtonian 
value" . 

k N 1 2 



2 12 24 

e. ^xi^^v^) 

ep, -.(- ~v-) ei-^C2 = - -%-/X''^ 

Iso _6_ -Ei-KEf-3B2)^^^ 

9+7J 3B2 



Iso -l + (-ci/e2)^/'' 

-Pfc C1+C2 
pole 4(3+T)) 



^Pk 36-35rj 
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TABLE 11: Taylor coefficients of the flux, phase, time and frequency. A'^ denotes the "Newtonian value" and 9 — 
[')(tiso-t)/(5m)]-i/8. In all cases the fe = coefficient is 1 and the k = 1 coefficient is zero. In certain cases the 2.5 
PN term involves log v or 6^ log 6 term rather than a or 6^ term. In those cases we conventionally include the 
logu dependence in the listed coefficient. Chirp parameters rt, k > 1, are related to the expansion parameters 
and (f>^ via Tk = (8(^fc — 5tj!)/3. We have given the simplified expressions for these in all cases, except k = 5 where no 
simplification occurs due to the presence of the log term in <j)^ . 

k N 2 3 4 5 

:p 32t)^ii^° 1247 _ 35t) ^ 44711 _,_ 9271i7 , 65 _ / 8191 _i_ 583j) \ 

5 336 12 9072 504 + 18 I 672 24 ,) ''^ 

fv 5m 743 , lln _ 32n_ 3058673 , 542977 _,_ 617i)^ _ ( 7729 _ 13 >i 

fe 256>7t>» 252 + 3 5 508032 ' 504 72 I 252 3 U 

Iv 1 3715 I 55t) _in._- 15293365 , 27145t7 , 3085»7^ / 38645 65?7 \ i / v \ 

T^fe 16>7t;=> 1008 ' 12 " 1016064 1008 ' 144 I 672 8 / visa J 

It _ 2 3715 I 5577 _3n 9275495 i 28487577 , 1855»7^ / 38645 6577 \ i / \ 

™fe rie^ 8064 96 4 14450688 + 258048 + 2048 1 21504 256 ) ^'^^ y g^^^ J 

pt 0^ 743 I II77 _37r 1855099 , 56975)7 , 37177^ _ ('JHSS. _ JJLnW 

fe Sirm 2688 + 32 10 14450688 + 258048 + 2048 21504 256 

^fe T^W l(W + 11^) -16^ 202 i (801 - 5t^) 



TABLE III: Fully maximised overlaps of the fiducial exact {X) waveform (effective-one-body signal [7]) with: (1) the 
standard time-domain post-Newtonian approximations of type tl and t3 (T*^, T*"^), given in Eq. ll^nj and Eq. (ITHll 
and Table HTl (2) the frequency-domain usual stationary phase approximations of type 1 and 2, (T*^, T^^) given in 
Eqs. 13.51 . 13.31 and 13.61 and Table UTI and (3) the time-domain P-approximants (P) - energy function as given 
by Eqs. THM and Eqs. ITBl . fl ux function in Eqs. 12.201 and 12.211 . and coefficients enlisted in Tables HI and HTl 
The overlaps, which are computed using the LIGO noise curve, are maximised not only over the time-lag and the 
initial phases of both the fiducial exact signal and the approximate template (by using two signal and two template 
waveforms, with phases equal to and 7r/2 [5]), but also over the two masses m\ and m2 of the approximate signal 
models. (The optimal masses are given below the overlaps.) The time-domain T*''-approximants are terminated 
when F = Q and the other signals are terminated when F(t) = /iso, the Iso frequency being determined consistently 
using E'jt^{v) — where Ea{v) is the corresponding approximate energy function. 



k 


(x,r,") (x,T«) 


{X,Tt) 


(x,rf) 




mi = m2 = 15M0 


4 


0.8881 0.9488 


0.8644 


0.8144 


0.8928 




(15.2,14.1) (16.3, 16.4) 


(14.7, 14.9) 


(16.0,16.1) 


(14.7, 15.1) 


5 


0.8794 0.8479 


0.7808 


0.8602 


0.8929 




(17.3, 16.4) (17.6, 17.9) 


(16.8, 16.7) 


(15.2,14.4) 


(15.4,14.3) 


mi = m2 = IOMq 


4 


0.9604 0.9298 


0.9581 


0.9109 


0.9616 




(10.1,9.6) (10.5, 10.3) 


(10.0, 9.7) 


(10.5, 10.6) 


(10.0, 10.2) 


5 


0.8814 0.8490 


0.8616 


0.9529 


0.9610 




(11.4, 10.6) (11.4, 11.7) ( 10.7, 11.0) 


(10.3, 9.7) 


(10.4, 9.7) 




mi = 


lOM0,m2 = 


1.4M0 




4 


0.9847 0.9673 


0.9835 


0.9721 


0.9937 




(1.27,11.1) (0.95,16.6) 


(1.27, 11.1) 


(0.96, 16.4) 


(1.35, 10.5) 


5 


0.9452 0.6811 


0.9394 


0.9922 


0.9941 




(0.82, 20.4) (1.11, 15.7) 


(0.82, 20.4) 


(1.34, 10.5) 


(1.37, 10.2) 


mi = 1712 ~ 1.4M0 


4 


0.8828 0.8538 


0.8830 


0.8503 


0.9719 




(1.40, 1.39) (1.42, 1.39) 


(1.40, 1.39) 


(1.44, 1.37) 


(1.47, 1.34) 


5 


0.8522 0.7643 


0.8522 


0.9994 


0.9945 




(1.46, 1.35) (1.43, 1.38) 


(1.46, 1.35) 


( 1.45, 1.35) (1.49, 1.32) 
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TABLE IV: One-sided noise power spectral densities (PSD) of initial interferometers, Sh{f)- For each detector the 
noise PSD is given in terms of a dimensionless frequency x = f /fo and rises steeply above a lower cutoff /s. 

Detector /,/Hz /p/Hz 10^^ x S'fe(3:)/Hz~^ 

[(3.4s)-3» + 34a;-i + 20(1 - a;" + x*/2)/{l + xy2)] 
9.00 [{4A9x)-^^ + 0.16a;"*-'^^ + 0.52 + 0.32x2] 
75.0 [x~^ + 13x-^ + 9(1 + x^)] 
3.24 [(6.23a;)-5 + 2x~^ + 1 + a;^] 
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400 
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10 100 

m (solar mass) 



FIG. 1: Signal-to-noise ratios in GEO, LIGO-I and VIRGO when using as Fourier-domain template the post- 
Newtonian model Eq. 13.61 (T-^^), truncated at the test-mass Fibo = 44OOM0/m Hz (in thin lines), compared to 
the optimal one obtained when the template coincides with the fiducial "exact" (effective one-body) signal (thick 
lines). As usual, we averaged over all the angles. The overlaps are maximised over the time lags, the phases, and 
the two individual masses mi and m2. The plots are jagged because we have computed the SNR numerically by 
first generating the fiducial "exact" waveform in the time-domain and then using its discrete Fourier transform in 
Eg. 15. 311 . The greater SNR achieved by effective one-body waveforms for higher masses, as compared to Fig 1 of 
DIS2, is due to the plunge phase present in these waveforms. Observe that the presence of the plunge phase in the 
latter significantly (up to a factor of 1.5) increases the SNR for masses m > 35Mq. Using the effective one-body 
templates will, therefore, enhance the search volume of the interferometric network by a factor of 3 or 4. 
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FIG. 2: The frequency evolution of the various approximate models is compared with the fiducial exact (10, 10)Mq 
model in the LIGO band. To indicate the effect on the overlap, we also plot the weighting function l/hn{f) for initial 
LIGO (not to scale), which is a measure of the detector's sensitivity. 



